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We analyze the first results on charged particle multi- 
plicity at RHIC in the conventional eikonal approach and in 
the framework of high density QCD. We extract the frac- 
tion F of the hadron multiplicity originating from "hard" 
(i.e. proportional to the number of binary collisions) pro- 
cesses; we find a strong growth of this fraction with energy: 
Ffyi = 56 GeV) ~ 22%, while F(y/s = 130 GeV) ~ 37%. 
This indicates a rapid increase in the density of the produced 
particles. We outline the predictions of high density QCD 
for the centrality, energy, and atomic number dependence of 
hadron production. Surprisingly, the predictions of the con- 
ventional eikonal approach and of high density QCD for cen- 
trality dependence of hadron multiplicity at y/s = 130 GeV 
appear very similar. 



I. INTRODUCTION 

Recently, the first results of the charged multiplicity 
measurement at RHIC were published by the PHOBOS 
Collaboration Q|. These data marked the beginning of 
collider era in the experiments with relativistic heavy 
ions, and already excited a controversy about the mecha- 
nism of multi-particle production at this new high-energy 
frontier |^-|(| . One of the focusing points of the discussion 
has been the role of saturation effects |||||. 

The purpose of the present paper is twofold. First, we 
analyze the experimental results in the framework of the 
eikonal Glauber approach which has been shown to de- 
scribe well the data on particle production in heavy ion 
collisions at SPS energies ( a detailed description of the 
formalism and further references can be found in Ref . ) . 
We then extract the fraction of the hadron multiplicity 
originating from the processes scaling with the number 
of binary collisions ( "hard processes" ) and make predic- 
tion for the centrality dependence. Second, we discuss 
whether or not the present data indicate the onset of sat- 
uration behavior, expected in high-density QCD [p| pd] . 
We briefly review the physical picture of parton satura- 
tion and make predictions for the centrality, energy, and 
atomic number dependence of particle production in nu- 
clear collisions; we then compare these to the predictions 
based on the conventional eikonal approach. 



It has been conjectured long time ago g] that at suf- 
ficiently high energy the particle production in nucleus- 
nucleus collisions will be dominated by hard processes. 
However, the gross features of particle production at 
CERN SPS energies were found to be approximately con- 
sistent (see, e.g., 12 1) with the scaling in the number 
of participants N vart , accommodated by the "wounded 
nucleon model" [flill Fl. There are deviations from this 
scaling; however they are not very large - the WA98 
Collaboration, for example, finds |l3|] the charged par- 
ticle pseudo-rapidity density at mid-rapidity scaling as 
~ Npart- Tne advent of RHIC has pushed the en- 
ergy of heavy ions to the new frontier, and the first re- 
sults on charged multiplicity jjj at two different energies 
(\fs = 56 GeV and y/s = 130 GeV) could already enable 
us to evaluate the relative importance of hard processes 
in heavy ion collisions at collider energies. 

The PHOBOS Collaboration presented a comparison 
of their results with the multiplicity measured in pp 
and pp collisions, which shows that the particle pseudo- 
rapidity density per participant increases by approxi- 
mately 70% near pseudo-rapidity r\ = in Au — Au colli- 
sions at \fs — 130 GeV. This is indicative of a significant 
contribution from hard processes^. We will now quantify 
this conclusion using the eikonal approach. 



II. THE EIKONAL APPROACH TO HADRON 
PRODUCTION 



Let us assume that the fraction x of the multiplicity 
n pp measured in pp collisions per unit of (pseudo)rapidity 
is due to "hard" processes, with the remaining fraction 
(1 — x) being from "soft" processes. The multiplicity 



1 (Multi) strange particle yields |15| show a remarkable dis- 
agreement with this trend, but here we will address only the 
total charged multiplicity. 

2 While it is true that hard processes scale with N co u, there 
may also be a component of the "soft" interaction with the 
same scaling, see e.g. [ fili|Jl7| ] ; the multiplicity analysis alone 
cannot distinguish between these two options, but they would 
lead to different predictions for the transverse momentum 
distributions. 
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in nuclear collisions will then also have two components: 
"soft", which we assume is proportional to the number 
of participants N part , and "hard", which is proportional 
to the number of binary collisions N co u: 



an 



x) n p 



+ x n pp (N coU ) 



(1) 



The same functional dependence has been used by Gyu- 
lassy and Wang in their analysis J^] based on the HIJING 
model. 

PHOBOS Collaboration quotes the following fit to the 
data on the pseudo-rapidity density of charged multiplic- 
ity in non-single diffractive pp interactions fl9| : n pp = 
2.5 - 0.25 ln(s) + 0.023 In 2 (s). When used at y/s = 130 
GeV, this formula gives n pp (y/s = 130 GeV) ~ 2.25; at 
y/s = 56 GeV one gets n pp (^ = 56 GeV) ~ 1.98. These 
are the values that will be used as an input to (^) in this 
paper. 

The shape of the multiplicity distribution at a given 
(pseudo)rapidity 77 can now be readily obtained in the 
usual way (see, e.g., JtJ): 



da 
dn 



d 2 bV(n;b) (1-P (&)), 



(2) 



where Po(b) is the probability of no interaction among 
the nuclei at a given impact parameter b: 



P {b) = (l-a NN T AB (b)) 



MS 



(3) 



where ctnn is the inelastic nucleon-nucleon cross sec- 
tion, and Tab (b) is the nuclear overlap function for the 
collision of nuclei with atomic numbers A and B; we 
have used the three-parameter Woods-Saxon nuclear 
density distributions ifTsfl . For ^/s = 56 GeV we have 
used <jnn — 37 ± 1 mb, while for yfs = 130 GeV 
@ n n = 41 ± 1 mb was chosen basing on the interpolation 
of existing pp and pp data |2(|. For y/s = 200 GeV, we 
use onn = 42 ± 1 mb. The correlation function V(n; b) 
is given by 



V(n;b) 



1 



\j2iran{b) 



exp 



(n — n(b)Y 
2an(b) 



(4) 



here n(b) is the mean multiplicity at a given impact pa- 
rameter b, evaluated analogously to Eq. ([!]) ; the formulae 
for the number of participants and the number of binary 
collisions can be found in Q . The parameter a describes 
the strength of fluctuations; we have chosen the value 
a = 0.6 which fits the data well. 

In Fig.l, we compare the resulting distributions for two 
values of x — and x = 0.09 to the PHOBOS experimen- 
tal distribution, measured in the interval 3 < \r)\ < 4.5. 
(The reason for presenting the result corresponding to 
this particular value of x will become clear later; we have 
performed calculations for a wide range of x.). One can 
see that both x — and x — 0.09 distributions describe 



the data quite well; the curve with x = (pure wounded 
nucleon model) fits the data somewhat better. This sug- 
gests that the multiplicity in the 3 < < 4.5 pseudo- 
rapidity interval is due to the "soft" processes. This is in 
accord both with the perturbative QCD picture, where 
the minijet production is peaked at mid-rapidity, and 
with the soft production approaches 16 17f] . On the other 



hand, since the shape of the distribution is described well 
in our approach and is not very sensitive to the rela- 
tive contributions of hard and soft components (we have 
checked that this is true up to x ~ 0.12), we can use this 
distribution as a good handle on centrality. 

The total Au-Au cross section computed in our ap- 
proach at y/s = 130 A GeV is a to t = 7.04 ± 0.05 barn. 
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FIG. 1. Charged multiplicity distribution at yfs = 130 A 
GeV; solid line (histogram) - PHOBOS result |lj ; dashed line 
- distribution corresponding to participant scaling (x = 0); 
dotted line - distribution corresponding to the 37% admixture 
of "hard" component in the multiplicity (x = 0.09); see text 
for details. 



Let us now establish the correspondence between cen- 
trality and the mean number of participants. The mean 
number of participants in a nucleus-nucleus collision with 
multiplicity n is defined as 



N p art{n) 



f d 2 bV(n;b)(l-P (b))N part (b) 
Jd^bV(n;b)(l-P (b)) ■ 



(5) 



where N part (b) is the mean number of participants in 
the Glauber approach. Analogous formula can be used 
to evaluate the mean number of collisions N co u(n) (for 
explicit expressions, we refer to |7j). 

If we now apply the 6% centrality cut (the cut applied 
by the PHOBOS Collaboration to extract the multiplic- 
ity in most "central" collisions in the pseudo-rapidity in- 
terval \r/\ < 1) to the distribution of Fig. 1 computed 
with x = 0, we can extract the mean number of partic- 
ipants and the mean number of collisions in the events 
with this centrality: 
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J dn/d 2 6 P(n;&)(l - P Q (b))N part (b) 

no 

J dnj d 2 bV(n;b)(l- P (b)) 



(6) 



where n > hq corresponds to the 6% centrality cut in 
nuclear collisions. 

The numbers deduced in this way are the following: 

(N part ) = 339 ± 2; 

(N coli ) = 1026 ±6; yfs = 130 GeV. (7) 



and 



(N part ) = 334 ±2; 
(N coU ) = 921 ± 6; 



Vs = 56 GeV. 



(8) 



The numbers of participants (Q, |^) are consistent with the 
ones quoted by the PHOBOS Collaboration: (N part ) — 
343 ± ^(stat)t 7 u (syst) for VI = 130 GeV and (N p 
330 ± A(stat)tl"(syst) for = 56 GeV. 
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For the sake of completeness, we give in Table 1 the re- 
sults of the calculations of mean number of participants 
and collisions for different fractions of the cross section 
and different energies. One can extract from this Table 
also the mean numbers corresponding to different cuts on 
centrality; for completeness, we list them in Table 2 to- 
gether with the corresponding average impact parameters 
and the mean densities of participants in the transverse 
plane for = 130 GeV. 

Table 1. The mean number of participants and bi- 
nary collisions corresponding to different centrality cuts 
in AuAu collisions at different energies, as computed in 
Glauber approach. 



centr. 




56 GeV 


v^ = 


130 GeV 


y/s = 200 GeV 


cut 


(Np ar t 


) (NcoU) 


(N part ) (N coU ) 


(Npart) 


(Ncou) 


0- 5 % 


342 


949 


344 


1053 


344 


1074 


0- 6 % 


336 


927 


339 


1028 


339 


1049 


- 10 % 


313 


847 


316 


937 


317 


958 


- 20 % 


266 


684 


268 


755 


270 


776 


- 30 % 


228 


563 


231 


622 


230 


634 


- 40 % 


196 


469 


198 


516 


199 


528 


- 50 % 


169 


393 


172 


434 


172 


444 


-100% 


92 


206 


93 


226 


93 


231 


10- 20 % 


218 


522 


221 


575 


222 


590 


20- 30 % 


150 


317 


153 


348 


153 


356 


30- 40 % 


101 


184 


102 


201 


102 


204 


40- 50 % 


63 


98 


64 


106 


64 


108 


50- 60 % 


37 


47 


38 


52 


38 


52 


60- 70 % 


19 


21 


20 


22 


20 


22 


70- 80 % 


8.7 


7.8 


9.4 


8.8 


9.5 


8.8 


80-100% 


2.5 


1.9 


2.7 


2.2 


2.8 


2.2 



We are now in a position to evaluate the "soft" and 
"hard" contributions to the observed multiplicity by us- 
ing Eq. (0) and the experimental values of dn/drj = 
555 ± U(stat) ± 35{syst) at ^/s = 130 GeV and dnjdj] = 
408 ± 12(stat) ± 30(syst) at y/s = 56 GeV. Using 
n pp = 2.25 as follows from the fit to the data JTgj ] , we get 
x = 0.09 ± 0.03 at s/s = 130 GeV; the use of the same 
fit at y/s — 56 GeV yields n pp — 1.98 and leads to the 
value x = 0.05 ± 0.03. These numbers are qualitatively 
consistent with the predictions based on the mini-jet pic- 
ture (see ]2l| ] for review), which indicate that the mini-jet 
contribution is increasing as a function of energy between 
y/s = 56 GeV and y/s = 130 GeV. 

Using Eq.((|), we find that at y/s = 56 GeV the frac- 
tion F ~ 22% of the produced particles result from hard 
processes, while at ^/s = 130 GeV this fraction increases 
to F ~ 37%; we define F as 



F = 



dn/drj 



(9) 



Table 2. The impact parameter dependence of the 
mean number of participants and binary collisions in 
AuAu system at = 130 GeV; we also list the average 
densities of participants in the transverse plane and the 
corresponding values of saturation scale Q s in the high 
density QCD approach computed according to Eq. (|l4|) . 



b 

(fm) 




NcoU 


Ppart 
(fm" 2 ) 


Ql 

(GeV 2 ) 


0. 


378.4 


1202.7 


3.06 


2.05 


1. 


372.4 


1173.6 


3.05 


2.04 


2. 


354.7 


1092.9 


3.01 


2.02 


3. 


327.0 


975.7 


2.95 


1.98 


4. 


292.2 


837.0 


2.86 


1.92 


5. 


253.2 


689.6 


2.75 


1.84 


6. 


212.3 


543.5 


2.60 


1.74 


7. 


171.5 


406.7 


2.41 


1.61 


8. 


132.4 


285.6 


2.17 


1.46 


9. 


96.5 


184.8 


1.89 


1.26 


10. 


65.1 


107.4 


1.54 


1.04 


11. 


39.5 


54.1 


1.15 


0.77 


12. 


20.5 


22.7 


0.72 




13. 


8.7 


7.8 


0.35 




14. 


2.9 


2.2 


0.12 




15. 


0.8 


0.5 


0.03 





While we have shown that the shape of the minimum 
bias distribution is not very sensitive to the relative mag- 
nitude of "soft" and "hard" contributions, the correlation 
between the forward energy and the multiplicity pro- 
duced at mid-rapidity is 0. To demonstrate this, we 
plot the forward energy Ep = (A — N par t/2)y/s/2 versus 
the multiplicity Eq. (Q) in Fig. |^. However if one nor- 
malizes the purely "soft" contribution to the measured 
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central multiplicity, the difference in the correlations be- 
comes less pronounced - see Fig. @. 

One can see that the correlation found in the case of the 
wounded nucleon model (x = 0) is significantly different 
from the correlation with x = 0.09, which corresponds to 
the fraction F ~ 37% of multiplicity produced in hard 
processes. Since at RHIC the zero degree calorimeter 
detects only the energy carried by neutral particles, the 
translation of the forward energy to the energy measured 
by the ZDC requires a dedicated study. 



700 




FIG. 2. Correlation between the charged multiplicity near 
r\ = and the forward energy; dash-dotted line - the cor- 
relation corresponding to the participant scaling; solid line - 
the correlation containing 37% admixture of "hard" compo- 
nent in the multiplicity; dashed line is the participant scaling 
correlation normalized to the measured central multiplicity. 
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FIG. 



3. Centrality dependence of the charged multiplic- 
ity per participant pair near 77 = at \fs = 130 A GeV; 
the curves represent the prediction based on the conventional 
eikonal approach, while the diamonds correspond to the high 
density QCD prediction (see text). The square indicates the 
pp multiplicity. 

The formula (Q) and the calculated numbers of partic- 
ipants and collisions given in Tables 1 and 2 allow us to 



compute also the centrality dependence of the multiplic- 
ity. With the parameters described above, formula (|l|) 
leads to 



2 dn 



1 



Np ar t dr) N, 



part 



(2.04 N part + 0.40 N coH ) 



(10) 



The result is shown in Fig. (^). It is qualitatively, but not 
quantitatively, consistent with predictions of the HIJING 
model presented in ||. 

A qualitatively different behavior has been predicted 
in the framework of the saturation approach of Ref . || , 
in which the multiplicity per participant decreases as a 
function of centrality. The approach || assumes the satu- 
ration of the produced partons, whereas the original sat- 
uration ideas concern the behavior of partons in 
the initial wave function of the nucleus. While the two 
approaches bear some similarities, they also have impor- 
tant differences. These differences, as we will see, result 
in qualitatively different predictions for the centrality de- 
pendence of particle multiplicity. It is therefore worth- 
while to examine the predictions of the (initial state) sat- 
uration approach. 



III. HIGH DENSITY QCD AND HADRON 
PRODUCTION IN NUCLEAR COLLISIONS 

Let us begin by sketching the basic ideas of satura- 
tion [^| |HJ in their most rudimentary form. Consider the 
wave function of the nucleus boosted to a large momen- 
tum. The nucleus is Lorentz-contracted, and partons 
"live" on a thin sheet in the transverse plane. Each par- 
ton occupies the transverse area tt/Q 2 determined, by un- 
certainty principle, by its transverse momentum Q, and 
can be probed with the cross section a ~ a s (Q 2 ) tt/Q 2 . 
On the other hand, the entire transverse area of the nu- 
cleus is Sa ~ ttR 2 a . Therefore, if the number of partons 
exceeds 



1 



(11) 



they will begin to overlap in the transverse plane and 
start interacting with each other, which prevents fur- 
ther growth of parton densities^]. This happens when 
the transverse momenta of the partons are on the order 
of 



Na 
r a 



(12) 



3 Closely related to this picture are the ideas of string perco- 
lation [E2lj23] in the transverse plane. 
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which is called the "saturation scale" . In the saturation 
regime, as is apparent from (11) and (13), the multiplicity 



1.23 ±0.20. 



of the produced partons should be proportional to 



Q\R\ ~N A ~A. 



(13) 



In the weak coupling regime, the density of partons 
becomes very large, which justifies the semi-classical 
McLerran-Venugopalan approach JToj - In the first ap- 
proximation, the multiplicity in this high density regime 
scales with the number of participants. There is, how- 
ever, an important logarithmic correction to this from 
the evolution of parton structure functions with the sat- 
uration scale Ql, which we discuss below. 

While our "derivation" has been very simplistic, the 
formulae (12 ^) can be reproduced by more sophisti- 
cated methods which also allow to reconstruct 



the coefficient of proportionality in (12): 



Q2 = 87T 2 N C ^ tr& 



a s (Qi) xG(x,Qi) 



(14) 



where iV c = 3 is the number of colors, xG(x,Ql) is the 
gluon structure function of the nucleon, and p par t is the 
density of participants in the transverse plane. We divide 
Ppart by 2 to get the density of those nucleons in a single 
nucleus which will participate in the collision at a given 
impact parameter. 

Let us estimate the saturation scale from ( |l4| ) for a 
central Au - Au collision at ^/s — 130 GeV. Eq.(p|) 
is an equation that can be solved by iterations; a self- 
consistent solution can be found at Q 2 S ~ 2 GeV 2 if we 
use xG(x,Q 2 ) ~ 2 j|| at x ~ 2Q s j^~s ~ 0.02, with 
a s (Q 2 )~0.6. 

As before, we will normalize the prediction to the ex- 
perimental data referring to the 6% of most central colli- 
sions. We use an explicit expression |llf| for the number 
of produced partons 



d 2 N 
cPbdrj 



N, 



1 JLga 



(15) 



where c is the "parton liberation" coefficient accounting 
for the transformation of virtual partons in the initial 
state to the on-shell partons in the final state. Integra- 
tion over the transverse coordinate and the use of ( |l4] ) 
yield simply 



dN 
dq 



= cN, 



part 



(16) 



If we assume that dN/drj ~ 3/2 dn exv jdr\^ take 
xG(x, Ql) ~ 2 at Ql ~ 2 GeV 2 and x ~ 2Q s /*/s~ ~ 0.02 
Hi, and use N part ~ 339 from Table 1 for the 6% 
centrality cut, the experimental number jlj] dn/drj = 
555 ± 12(stat) ± 35(syst) translates into the following 
value of the "parton liberation" coefficient: 



(17) 



This number appears to be close to unity, as expected by 
Mueller ]Tl[ | , which implies a very direct correspondence 
between the number of the partons in the initial and fi- 
nal states. Moreover, this may imply that the number 
of particles is conserved through the parton-to-hadron 
transformation - a miraculous fact first noted in the con- 
text of the "local parton-hadron duality" hypothesis [E5| . 

The value (O) can be compared to the recent lat- 
tice calculationj^6j by Krasnitz and Venugopalan, which 
yields c = 1.29 ±0.09. Very recently, an analytical calcu- 
lation for c has been presented by Kovchegov 27 1, with 
the result c = 2 In 2 ~ 1.39. 

To compute the centrality dependence, we still need to 
know the evolution of the gluon structure function with 
the density of partons, which is proportional to the mean 
density of participants in the transverse plane. We will 
assume that this evolution is governed by the DGLAP 
equation 28 2^] , and take 



In 



A 2 

QCD 



(18) 



The dependence (18) emerges when the radiation of glu- 
ons is treated classically and so is consistent with Eq. 
(|l4|). Eqs. (|lj) and can now be used to evaluate 
the centrality dependence; with the parameters described 
above we get 



N, 



part 



dn 
dr) 



0.82 In 



Ql 



A 2 

QCD 



(19) 



where we take Aqcd — 200 MeV and use the values of 
Q s listed in Table 2. 

We present the results of our calculations in the satura- 
tion scenario in Fig. [| The similarity of the predictions 
based on the conventional eikonal approach and on the 
high density QCD is striking in spite of the a priori to- 
tally different functional dependences (|lfj| ) and (|l9|). 



IV. SUMMARY AND DISCUSSION 

To summarize, our analysis of the first RHIC results 
M implies that the role of hard processes in particle pro- 
duction rapidly increases at collider energies, bringing us 
to the regime of high parton densities. The experimen- 
tal study of centrality and atomic number dependence of 
particle production at different energies is needed to ex- 
tract the crucial information about the behavior of QCD 
at these conditions. 

We have made predictions for the centrality depen- 
dence of the particle multiplicity both in conventional 
eikonal approach and in the framework of high density 
QCD. The latter prediction is different from the ones 
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available in the literature |3|,g] and shows that the multi- 
plicity per participant increases as a function of central- 
ity, reflecting the evolution of parton densities with the 
increasing saturation scale. Surprisingly, both conven- 
tional and high density QCD approaches lead to almost 
identical dependence on centrality. This makes it dif- 
ficult to distinguish between the two approaches using 
only the data on centrality dependence of multiplicity at 
one fixed energy, and one has to rely on the analysis of 
the transverse momentum distributions. We leave this 
for the future. On the other hand, the numerical simi- 
larity of the predictions of eikonal and high density QCD 
approaches means that the apparent success of the con- 
ventional approach in describing the previously available 
data might mask a different physics. 

The approach based on the saturation allows us to ex- 
tract the initial energy density of partons achieved in 
AuAu collisions at RHIC. Indeed, in this approach the 
formation time of partons is tq ~ l/Q Sl and the trans- 
verse momenta of partons are about kt — Q s . We thus 
can use the Bjorken formula Q and the set of parame- 
ters deduced above to estimate 



< k t > d 2 N 



To 



cPbdr) 



d 2 N 
d?bdr\ 



18 GeV/fm 3 . (20) 



This value is well above the energy density needed to 
induce the QCD phase transition according to the lattice 
calculations. 

High density QCD also provides very definite predic- 
tions for the energy and atomic number dependence of 
hadron multiplicity - according to (16), the multiplicity 
simply scales with the number of participants and the 
gluon structure function of the nucleon, which one has 
to evolve to the saturation scale. The knowledge of the 
extracted value ([l7|) furnishes the prediction. 

We are fully aware of the perils involved in applying 
the weak coupling methods at the scales on the order 
of Ql ~ 1 GeV 2 . There are reasons to believe (see, e.g., 
JIi| ) that the gluon distributions at these scales may have 
a strong non-perturbative component, which would in- 
fluence the results. The clarification of related problems, 
as well as the problem of subsequent evolution of the 
produced partons, await further work. 

We hope that the forthcoming experimental results 
and further progress in theory will eventually allow to 
uncover the dynamics of QCD at the high parton den- 
sity, strong color field frontier. 
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A.H.Mueller, H. Satz and X.-N. Wang for helpful conver- 
sations on the issues related to this work. The research 
of D.K. is supported by the U.S. Department of Energy 
under Contract No. DE-AC02-98CH10886. 
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